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Abstract 

A global analysis of the full amplitudes for e + e~ — > VP (Vector and Pseudoscalar) channels at yfs 
=3.773 GeV and 3.670 GeV, which were measured by the CLEO-c Collaboration, shows that those 
measurements are essentially nontrivial for searching for the ■0(3770) non-DL> decays. Unlike the 
nearly negative verdict on the 0(3770) strong decays to the VP channels in the original analysis of 
the CLEO-c data, there exist some unusual solutions that predict the remarkable strength of SU(3) 
symmetry VP decay of 0(3770) resonance, which give some clue to understand the mechanism of 
■0(3770) non-DD decays and to reexplain the well-known p — tt puzzle in the J/0 and 0(3686) 
decays. 
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I. MOTIVATION 



There is a long-standing puzzle in understanding the exist measurements for ^(3770) 
and DD production cross sections at the peak of ?/>(3770) production in e + e~ annihilation 
[l|, 2)]. Potential Models predict that ?/>(3770) decays into DD with branching fraction of 
~ 100%. Recently careful investigation shows that the branching fraction of ^(3770) non- 
DD decay would be up to more than 10% 0, 4]. It is very interesting to know what are the 
exclusive non— DD final states of ?/>(3770) decays. Except the electromagnetic transitions 
and hadronic transitions of ?/>(3770) to lower charmonium states, are there indeed other 
significant exclusive non-charmed decay modes from ■0(3770) decays? 

In the charmonium decays, there is another long-standing puzzle in understanding the 
p7r decays of J/ip and ?/>(3686). The partial widths of pre channel and other VP channels in 
the ^(3686) decays are unexpectedly lower than those in J/ip decays. This is so called "p7r" 
puzzle. Are the J/ip decay rates enhanced by some unknown mechanism or the ^(3686) 
decay rates are suppressed abnormally? To investigate the possible source of this puzzle, it 
is also important to measure the ^(3770) VP decay amplitude. 

Recently, the BES Collaboration |5| observed a large production cross section for e + e~ — > 
K*(892)°K°+c.c. 



a(e + e- K*°K + c.c.) = (15.0 ± 4.6 ± 3.3) pb, 

at center-of-mass energy of ^=3.773 GeV and found that the K* ± (892)K T production is 
suppressed. Taking into account the possible interference between the strong decay ampli- 
tude and the continuum production amplitude at a/s=3.773 GeV, the BES Collaboration 
set an upper limit on the strong decay partial width for ^(3770) — > K*(892)K+c.c. to be 

T(^(3770) ^ K*(892)K + c.c.) < 29.0 keV 

at 90% confidence level. 

The CLEO-c Collaboration made more careful studies of twelve exclusive VP decay chan- 
nels for ^(3770) — > pn, K*(892)K + c.c, cure , pi], pi]', ui], ui]', (pi], (pi]' and 07r° reported 
in Ref. [6J. The CLEO-c Collaboration measured the cross sections for all of the channels 
at the energies y / i=3.773 GeV and a/s=3.670 GeV. The CLEO-c results show that the 
measured cross sections at -^5=3. 773 GeV are almost equal to or even less than the ones 
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measured at y / s=3.670 GeV, which mean that the net cross sections for the ^(3770) decays 
are consistent with zero except only for the channel ^(3770) — > (jnj. The negative results 
about the ^(3770) strong VP decays led people ignoring the important strong decay com- 
ponent existing in ^(3770) and only focusing their attention on the form factors of those 
channels as well as the isospin violation in electromagnetic interaction 

In this paper, we develop a model to account for both the amplitudes of electromagnetic 
(E-M) production and t/>(3770) strong decay in the process of e + e~ — > VP. By analyzing 
the cross sections for the exclusive VP channels, which were measured by the CLEO-c 
Collaboration, we extract out the branching fractions for ^(3770) decay to these VP final 
states. 



II. THE MODEL AND THE FORMULAE 



In the ?/>(3686) decay sector, because of the smallness of the strong VP decay coupling, 
the E-M decay component as well as the continuum (E-M) component of the VP channel 
would be no longer the small amounts comparing with those of strong decay. People have 
to deal with the two components properly {t|. At the resonance peak, the production 
amplitude consists of two parts, one is the decay amplitude of charmonium resonance and 
the other is continuum E-M production amplitude. In the resonance decay part, there 
are two components as well. They are the E-M decay amplitudes and the strong decay 
amplitude. Totally, there are three components involved in the e + e~ annihilation process at 
the resonance peak, which are the strong decay component, the E-M decay component and 
the continuum production component. 

Unlike the VP decays of J/ip and ^(3686), the E-M decay amplitudes of ^(3770) can 
be neglected due to the little tiny dileptonic decay branching fraction. There are only 
the strong decay amplitude and the continuum production amplitude in the e + e~ collision 
production at -y/i=3.773 GeV. Typically, according to the conventional point of view, the 
partial widths of the ^(3770) VP decay channels could be up to keV order of magnitude, 
like their cousins in J/ip decays. However, due to the large width of ^(3770) the decay 
amplitudes of those channels can not get large amplification as the ones at the narrow 
resonance states. Associated with the measurements of the form factors of channel utt at 
the energies of -^(3686) and -^(3770) resonance vicinities js, [§], [9] the decay process with 
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only a few keV partial width of the rather wide resonance is really hard to be measured if 
one does not consider the interference between the amplitudes of the strong decay and the 
continuum production. As the measurements by CLEO-c [6], both the evident yield excess 
of channel <f)7] and the rather large yield deficit of the pir channel at the resonance peak 
hint that there must be rather complex interference between the two kinds of amplitudes 
acting globally on the VP channels. Some destructive interference just shows up at pn 
channel in the "deficit" way. And more complex interferences cause the <fir) yields enhanced 
at resonance peak. In fact, in such complicate interference case the decay contributions may 
easily be covered up by the continuum contribution and the interference contributions. If 
one completely neglects the buried decay contribution, the single E-M amplitude assumption 
would not describe the measured cross sections well. In practice, it is dangerous to measure 
the branching fractions for the ?/>(3770) non-DD decays by simply considering the net yields 
for the channels observed at the peak of t/>(3770) over that at the nearby off resonance region. 
In this analysis we introduce the strong decay amplitudes in the analysis formalism to see 
how the strong decay affects the VP production at the -^(3770) resonance peak. 

We describe the global decay of ?/>(3770) and the continuum production process still based 
on the flavor SU(3) invariant model, which was developed thirty years ago 
the strange quark mass correction in both of the strong coupling and the E-M coupling, 
the wave function nonet symmetry breaking and the double Okubo-Zweig-Iizuka (DOZI) 
suppression effects are all taken into account. As for the continuum production at the two 
energy points v^=3.773 GeV and 3.670 GeV, except the coherent strong decay amplitudes 
from the J/tp and ?/>(3686) tails, which can safely be neglected from the calculations, there 
is only the continuum E-M amplitudes itself. In addition, we guess that the incoherent 
component contributions which are mainly from the initial state radiative (ISR) return to 



10] , In this model 



J/ip and ?/>(3686) resonances have efficiently been rejected in the work reported in Ref. js] 
and can be neglected in our analysis too. 



Following the convention given in Ref. Uj, we define that g represents the VP strong de- 



cay amplitude in the flavor SU(3) symmetry limit; g s represents the strong decay amplitude 
from s quark, 

s g = 1 - (&/|g|)/2, (1) 
characterizes the 577(3) mass violation, which is as the same as the parameter "s" given in 



the Tab. VIII of Ref. 



11] ; 6p represents the 77 — 77' mixing angle; the product r ■ g represents 



the amplitude correction of the SU(3) nonet symmetry violation with the factor (1 — sp) 
for a strange pseudoscalar production and with the factor (1 — sy) for a strange vector 



production. If sy = sp=0, (exactly sy 
correction. Unlike the case in Ref 



sp=0), r • g measures the pure DOZI amplitude 
111 ], because the E-M amplitude is no longer small 
comparing with the strong amplitude, we have to consider both the isoscalar and isovector 
components of the E-M amplitude. We define the E-M amplitude in the form of SU(3) octet 
matrix representation as 

E = ei • I 3 + e • Y 

in which e and ei are the isoscalar and isovector components, respectively, and I3 and Y 
are, respectively, the isospin third component and the hypercharge matrices in flavor SU(3) 
octet space. We define 9q as the phase of e relative to g, 8\ as the phase shift difference of 
ei to e and a factor (1/2 — s e ) as the correction for strange quark coupling to E-M isoscalar 
part e . We assume that the couplings e's and their phases do not change in the all VP 
channels, and their moduli at the two different energy points y / 5=3.773 GeV and a/s=3.671 
GeV only change with a 1/s 3 dependence. If e x = e , we return to common definition as 



Refs. 
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In! did. 

For the channels "c/i", (ch= pV,rV + c.c, etc) at energy a/5, M^ s ^ denotes the 
resonance decay amplitude and M^ m ^ denotes the continuum production amplitude. The 
total production amplitudes are then written as 

M$ = M^ s + M c t^, (2) 

in which 

M ch — n 



and 

1\/Tch A/T ch f 

iw ctm,3670 ~~ lvl ctm,3770 ' Jd 

at the energy of 1/5=3.670 GeV, where fd = 3.773 3 /3.670 3 is the scaling factor for the 1/s 3 
energy dependence of the cross section. The amplitudes M^ m377Q and Mfl s 

,3770 defined at 

^=3. 773 GeV for all of the channels are listed in Tab. HI In the table, X v = cos(54.736° + 
Op), Y v = sin(54.736° + fl P ), Xj = -sin(54.736° + ^ P ) and Y v > = cos(54.736° + ^ P ), which are 



the same as those given in Ref. Uj. Those amplitudes completely control the correlations 
among the VP channel productions. If any significant decay amplitude M^ s 3773 given in 
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TABLE I: The Amplitudes for VP production in e + e~ annihilation at yfs = 3.773 GeV. The 
coupling gK defined in channels K*K + c.c. can be considered as a free parameter if one of g and 
g s is fixed. 



Channel(c/i) 


M res,3770 




-^ctm.,3770 




g 




eo 


ujrj 


gX v + V2rg[V2X v + (1 - s P 


)^] 




4>V 


g s K, + rg(l- S v)[\ / 2A„ + (l- 


sp)^] 


-2e (l - s e )F„ 


UJT]' 


gXrf + V2rg{V2X v , + (1 - s P 






(f>rj' 


g s Y v , + rg(l - s Y )[y/2X vl + (1 - 




-2e (l - s e )Y n > 









3ei 












p°r] 







3eiX v 


p V ' 







3e x Y v 


K*°K° + c.c. 


SK = (g + g s )/2 




-eo(l/2 - s e ) - 3/2ei 




gx = (g + gs)/2 




-e (l/2 - s e ) + 3/2ei 



Tab. |T]has been measured to be non-zero, which means that the measured couplings g, g s 
etc. are non-zero, this indicates that ?/>(3770) has a significant branching fraction for decay 
to the non-charmed channel u ch" . 

We start the analysis from the observed numbers of the events for the VP channels, their 
errors and their corresponding detection efficiencies at the energies of ->/i=3.773 GeV and 



3.670 GeV, which were published by CLEO-c Collaboration [6]. The numbers and errors of 
the events are obtained in both the signal windows and the side bands. Tab. HT1 shows those 
numbers and detection efficiencies which are given in Ref. {(j). Taking the numbers of the 
events from Tab. HI1 we obtain the numbers, 

N o^ch = {N g _ N Vs u (3) 

of the observed events at y/s for the channel u ch v . Usually the determinations of the detec- 
tion efficiencies and the ISR corrections are all energy dependent and relate to the production 
line shapes for those channels. We guess that the determinations of the detection efficiencies 
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TABLE II: The numbers of the events observed by CLEO-c, the subscripts sw and sb indicate the 
signal window and the side band window, respectively; the upper script 3.67 and 3.77 indicate the 
cm. energies. 



Channel(c/i) 


N 3.67 
x " sw 


^3.67 


AT3-77 

- 1 ' sw 


A^ 7? 


<%) 


pn 


43 


5.4 


314 


44.8 


26.3 


A 


21 


3.4 


130 


33.0 


32.5 




22 


2.0 


184 


11.8 


23.1 


ooir 


54 


6.2 


696 


39.2 


19.0 




1 


1.6 


2 


40 


16.5 


p°rj 


36 


3.1 


508 


31.0 


19.6 


LOT] 


4 


0.0 


15 


6.0 


9.9 




5 


1.0 


132 


15.9 


11.0 


pV 


1 


0.0 


27 


0.9 


2.9 


LOT]' 





0.0 


2 


0.0 


1.5 


H 





0.0 


9 


2.0 


1.2 


K*°K° + c.c. 


38 


0.4 


501 


18.1 


8.8 




4 


1.0 


36 


32.4 


16.0 



given in Ref. js] were done under the assumption of that the continuum cross section line 
shape is in the 1/s 3 energy dependence and the energy cut is at vs 7 > J/ip mass, where 
is the center of mass energy of the ISR return system. It should be stressed that, for 
the ISR and FSR (Final State Radiative) corrections in the continuum processes, our cal- 
culation gives r] ct m=l.l9 at -^=3773 GeV and ?7 rfm =l.ll at ^=3670 GeV, while Ref. js] 
gives rjct m =l/l. 20=0. 833 at both of the two energy points. As for the resonance decay, the 
ISR correction is quite different from the one for the continuum process. In our calculation, 
the ISR correction factor for the resonance is rj res =0.824 including the FSR correction at \/i 
=3.773 GeV for all channels. In our ISR correction calculations, the v 3 phase space depen- 



dences have been taken into account, where v = ^[l — {my + mp) 2 /s][(l — (my — mp) 2 /s] 
is the velocity of the vector daughter in the CM decay system. However, in the calcula- 
tion of ^ctm's, a mean production threshold of the channels has been set to serve as the 
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common threshold for all of those channels. For this reason, the channel dependences of 
the corrections r^m's are ignored in this analysis. For the calculation of the contribution 
of the interference terms among the amplitudes describing different processes, we have to 
know their own detection efficiencies and ISR corrections. In the analysis we simply take 
the geometric average of the related coefficients as the effective ones. 

Using the amplitudes of VP channels in Eq.(J2]) and Tab. [J the efficiencies for those 
channels, the ISR and FSR correction as given above, the v 3 phase space dependences and 
the luminosities accumulated at the two collision energy points of yfs =3.773 GeV and 
3.670 GeV, we can calculate the expected numbers N^q cH and N^fy of the event yields 
for those channels under the simplification assumption concerning the interference terms 
mentioned above. For channel 11 ch" 

AT exp,ch _ j 3 ch ch i M ch |2 

iv 3670 ~~ -^3670 ' y 3670 ' t ctm ' '/cim,3670 ' l Jvi rfm,3670l 

and 

'> ch — T „,3 



^3773° — -^3773 ' "^3773 X 



. / c ch . rr,ch . njch _i_ ch ^ch ]\/f 

V t res 'Ires lvl res, 3773 ~T~ V ctm ' ctm < 



rch 

L ctm,3773 



Comparing the numbers N°^' defined in Eq.(J3]) with the expected one N^P> ch ^ we g e t 
the equation set 

jyexp,ch jqobs,ch 

\/~s ^fs W 

in which ch = pn, K*(892)K+c.c, uti , pi], pi]', u)t),u)t)' ,<jyq, <prj' and (f)Tr°, y/s=3. 670 and 
3.773 GeV. Because of the zero observation and zero expectation for the 07r° channel, listed 
in Tab. Hi] and Tab. HJ we can get rid of this channel in our analysis. So we only focus our 
attention on the rest eleven channels. 

In the maximum likelihood fit, leaving the parameters g, g s , e , e 1; r, s e , sy, sp, 6p Si 
and cos#o free, we can solve the Eq.(jH) by maximizing the probability function 

2N ch 

Prob= \{ P ch (N°£ ch ,N e J?> ch ) (5) 

ch,y/s 

where P c f l (N'^' ch , N^? ,ch ) is the probability of finding ]\[°^f> ch events with the assumed mean 
number j\f^?> ch a t the energy a/s, and get the solution of Eq. (j5J) with most probable values 
of the parameters which control the VP channel production at i/i=3.773 and 3.670 GeV. 
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From the solutions of g, r, eo and/or ei etc, one can get the Born cross sections for the 
VP channel production at ^=3.773 or 3.670 GeV. For example, the Born cross section for 
decay channel ?/>(3770) — > "ch" at y/s =3.773 GeV can be written as 

ch 3 I \jjch 1 2 ia\ 

a res — y 3773 ' \ Ivl res,3773\ ■ \ U ) 

The total Born cross section for e + e~ annihilation to the channel "c/i" at -y/i=3.773 GeV is 
given by 

°T = ^3773 ' l^res,3773 + ^rfm,3773| 2 ' CO 

While the Born cross section for the channel u ch" in continuum production at y/s is then 
given by 



a 



" -^-IKLJ 2 - (8) 



ctm,y/s \/s I ctm,y/s 



III. THE RESULTS 



According to the different coupling configurations of the amplitudes for the VP channels, 
shown in Eq.(T5]) and Tab. [U the twelve VP channels can be divided into three sub-sets. The 
first one consists of the channels without E-M isovector components, such as pn, cur], (pi], urj' 
and <fir)'. The second one consists of the E-M isovector component only, which is the pure 
E-M channels um°, prj, and pi]'. The third subset includes only the channels K*K + c.c, 
for which the amplitudes involve all of the two E-M coupling parts and the strong coupling 
component as well as their interferences in the VP production. Since there is no common 
coupling parameter in the first two sub-sets despite of the pseudoscalar mixing angle dp, 
one can simply try to solve Eq.flU) separately in the two sub-sets at first. 



A. The channels without isovector E-M amplitude 

We start the analysis from the first set in which the channels are without the isovector E- 
M amplitude contribution. From the numbers of events, N°J' C , observed at the two energy 
points of y / i=3.671 GeV and 3.773 GeV for the six channels p°7r°, p ± n T , ur], (fir], ur/' and 
cfir}', we solve the Eq.Q. Leaving all of the related parameters, such as |g|, g s , |e |, cos9 , 
s e , Op, r, sy and sp free, the fitting yields 
|g| = 2.752 ±0.291, 
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g s = 1.279 ±0.706, 
cosflo = -0.935 ±0.141, 
|e | = 1.495 ±0.126, 
s e = 0.264 ±0.247, 
r = -0.236 ±0.070, 
6 P = (-17.83 ± 12.12)°, 
s v = 0.092 ± 0.941 
and 

Sp = -0.091 ±0.587. 

n 

Indeed, as hinted by the measurements from the CLEO-c [6|] and guessed above, the solution 
of cos#o = —0.935 ± 0.141 shows that g and e are almost opposite in the VP production 
at ^(3770) resonance peak. Associated with earlier measurements in the J/i/j production 
and decays, in which the E-M decay amplitudes and the strong decay amplitude are more 



likely with the phase difference of (j) e res — <Pg ~ ±90° 
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131 ]. and the phase difference 



between the continuum amplitude and the one of resonance E-M decay at the resonance 

n n 

peak is 4> e ,ctm — 4>e,res = ±90° too [U, 115], the measurement of the phase difference 6 = 
4>e,ctm ~ 4>g = (I6OI17) between e and g here is reasonable. The 90° phase difference 
between the E-M decay amplitudes and the strong decay amplitude shown in J/ip decays 



12 



131 ] is mainly due to the exist of an original 90° phase difference from the short distance 
force range in the on-shell three gluon annihilation of 1 quarkonium states. This argument 
should be kept in the -?/>(3770) decays. Out of the hard quark-gluon interaction level the 



strong phase shift of long distance final state interaction should essentially be small 16]. 
The resolutions of sy = 0.092 ± 0.941 and sp = —0.091 ± 0.587, which are consistent with 
zero, seem to mean the small nonet symmetry breaking between the singlet and octet wave 
functions. So we can assume that the nonet symmetry maintains for the wave functions, 
which means that the nonet symmetry violation is only in the DOZI coupling. If we fix sy= 
sp =0 in the fit, we get the parameters |g|, g s , |e |, cos9 Q , s e , dp and r, which are listed in 
the second column ("Solution 1") of Tab. [HH where 6 P = (-18.73^)° indicates 

X v = Y v , = 0.809l°;°6l and = -X^ = -0.588l°$8 
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TABLE III: The fitted parameters. 







Solution 1 


Solution 2 


Solution 3 


Solution 4 






the channels 


the channels 


all channels 


all channels 






without isovector 


with isovector 


with 


without 






E-M component 


E-M component counter terms counter terms 


Is] 




2 64 +cul 


2.64(fixed) 


9 ^7+0.38 
^•°'-0.32 


2 43+ - 31 

z -^°-0.29 






1 no+0.67 
l.zo_ .74 


_i 70+0.42 

L - ' y -0.36 


1.18_ .78 


1 90+O.52 
1 -^°-0.49 


S 9 








0.28 ±0.16 


0.75 ± 0.42 


p. add 

OS 








_1 r;i+0.43 
1 - ol -0.40 


O.O(fixed) 


cos#o 




-0-91±g : g 


-0.91(fixed) 


-0.91 ±0.07 


-0.861852 


l e o| 






1.49(fixed) 


1 zLQ+ 012 


1.44!^ 


Se 




n 25+ - 19 

u.zo_ 22 


-0.35 ±0.11 


20 +0 - 20 

U.ZU_ 023 


-0 29+ - 13 


a» 








_n c;c+2.38 

u.oo_ 22 


O.O(fixed) 


l e i| 






1 990+0.016 

i.zzo_ 018 


1 O9R+0.017 

±.zzo_ 019 


1.22 ±0.02 


61 






(7.39i 5 7 ;^)° 




(7.62±^|)° 


Op 




(-18.711;?)° 


(-23.2 ±2.2)° 


(-22.51^)° 


(-24.ilii)° 


r 




-0.24 ±0.06 




-0.26 ±0.06 


n 09+0.09 

u - oz -0.08 


s v 




O.O(fixed) 




O.O(fixed) 


— f) 42+ - 36 
u - 4Z -0.59 


s P 




O.O(fixed) 




O.O(fixed) 


0.45 ± 0.70 


X 2 /ndof 




4.63/5=0.93 


3.98/5=0.80 


9.26/11=0.84 13.2/11=1.20 



The fit gives x 2 / n dof = 4.63/5 = 0.93, which is also listed in Tab. IIHI With the g and g s , 
we obtain the 577(3) mass correction 

s - 2fi8 +0 - 140 

Except for the three coupling strengths, |eo|, |g| and |g s |, which have their own dynamics 
in the higher energy position, the other relative correction parameters, s g , s e , r and the 
mixing angle #p obtained from the fit are all reasonable comparing with those obtained from 
J /if) decays measured by Mark-Ill and DM2 [171 ] . However, from Eqs.flS]) and (Q, we find 
that the strong decay coupling |g| gives quite large cross section, branching fraction and the 
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partial width for 0(3770) — > pn decay, which are 

<C = (18.21^) pb 



B p7T = (1.83i£g?) x 10" 
PPTT = 49.7+16-9 keV _ 



(9) 



The partial width is almost in two order of magnitude higher than that of the conventional 
typical partial width of the J/ip VP decays. The latter one is at the order of 1 keV. If 
we assume that the fraction of the width of -0(3770) — > pn to the width of -0(3770) — > 
light hadrons is roughly as the same as the one in the J/-0 decays, the huge partial width 
of -0(3770) — > pn would predict that about 10% of -0(3770) decays to noia-DD final states. 
This predicted branching fraction for -0(3770) — > non—DD is almost as the same as the one 
measured by BES-II Collaboration 0, 0] . The large cross section a^ s as shown in Eq. (j9]) is 
a factor of more than 3 of the total pn production cross section, 

< = 5.351^ pb, (10) 

obtained from Eq.([7]). This cross section is consistent with the measurement reported in 
Ref. [6|, (see Tab. II Vp . In the Tab. IIVI we list all of the production cross sections and the 
branching fractions for ■0(3770) — > VP predicted in this work, and also listed the production 



cross sections of the VP channels given in Ref. js] as the comparison. As for the pir 
production cross section at y/s = 3.670 GeV in this measurement, we get 

<3670 = 6-341^ pb, (11) 



which is indeed higher than the one given in Eq.f jTOj) at resonance peak. Because of the 
different determinations of the ISR corrections in Ref. [6] and in our work, the decrease of the 
Born cross section of pn channel at 0(3770) resonance peak is not so large as that obtained 
by CLEO-c Collaboration [(J (see also Tab. IIVI) . Owing to the large cancellation between 
the two amplitudes g and e , the large cross section a^ s for 0(3770) — > pn disappeared 
without the global amplitude analysis. 

B. The pure E-M channels um°, p°rj and p°7]' plus the channel K*K-\-c.c 



We consider together the last two sub-sets of the channels in which the amplitudes of 
E-M production contain isovector component. The isoscalar E-M component and the strong 
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decay coupling only serve in the channel K* K+c.c, and the strong coupling here is with a 
combined form = (g + gs)/2. Inserting the ten numbers of the events observed at the 
two energy points of 3.670 GeV and 3.773 GeV for the rest five channels into Eq.flU), leaving 
ei, Si, s e , 6p and gx free (or instead of gx, leaving g s free but fixing g at some reasonable 
value such as |g|=2.635 obtained from last solution independently), and fixing |eo|=1.494 
and cos#o=-0.907 obtained also from last solution in assumption of that there is no more 
correction added to the couplings g and |eo| measured in last solution in the channel K*K, 
we fit the numbers of events observed in the five channels and obtain the solution of the free 
parameters. The results are listed in the third column ("Solution 2") of Tab. IIHI The fit 
gives x 2 /n d of = 3.98/5 = 0.80. 

From above two solutions, we see that the measured pseudoscalar mixing angles 9-p in 
the two independent measurements are consistent with each other. The isovector E-M com- 
ponent, ei, is really split from the isoscalar one, eo, with almost 2a deviation in magnitude 
and with a non-zero phase shift difference. However, we note that the s quark strong decay 
coupling g s and its E-M coupling correction s e in the solution 2 are both with the negative 
values, while they are supposed to be positive in the conventional 577(3) invariant model 
with simple static mass corrections. The minus s e is formally very likely to be some anoma- 
lous "magnetic moments" term added to the s quark E-M coupling. As for the minus g s , 
which is obviously irrelative to the three pure E-M channels, it seems that the s quark strong 
coupling undergo almost 180° phase shift from the conventional S77(3) strong interaction 
wave function. The odd behavior of the minus g s as well as the E-M "anomalous magnetic 
moments" term for the s quark indicate that there might be some other dynamic sources 
or more complicated interaction correction contributing to the production of K* K+c.c. in 
e + e~ annihilation. 

Owing to the cancellation of g and the opposite g s (small g#), the decay cross section, 
branching fraction and partial width 



<J K * K 

res 


= (0.558™) pb ' 






B K*K 


= (0.561°;™) x 10- 4 


> 


(12) 




= {1.52*%) keV, t 







for ^(3770) — > K*K+c.c in this solution are quite small comparing with that of ^(3770) — > 
pn measured in last solution given in Eq.flH]). These results are also listed in Tab. [TV] 
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The remarkable increase of the E-M coupling eo by negative correction s e and the large 
interferences between the two E-M amplitudes and the strong amplitudes result a serious 
asymmetry between the total production cross sections of channels K*°K° +c.c. and K*^K^ . 
For example, at y/s = 3.773 GeV the cross sections 



a 



K*°K +c.c. 
T 



19.30j*gj 



pb 



a 



< 0.50 (pb, at 90% confidence level), 



(13) 



are consistent with the observed values reported in Refs. and fl. Using the parameters 
|ei| and Op, and from Eq.flBJ) we get the production cross sections of the three pure E-M 
channels, 



a 



ctrn 



'11.74; 



-0.31a 
-0.34 J 



pb 



a 



pi) 



ctm 



(14) 



(8.00±0.22) pb 
a± = (2.57±0.28) pb.J 
These are also listed in Tab. IIVI The E-M production cross sections and the total production 
cross sections obtained in this subsection are all systematically lower than those measured 
by CLEO-c Collaboration [6] by about 30%, (see Tab. HV|) . Those differences are also due 
to the different determinations of ISR corrections in the two works as mentioned above. 



C. The global fit including all of the measured VP channels 

We can introduce two additional effective counter terms, g^ dd and a M , to compensate the 
odd behavior appeared in the channel K*(892)K+c.c. for both the strong coupling and 
E-M coupling of s quark. We can simply assume that g^ dd and g, a M and e are collinear, 
respectively. In this case, we define the amplitudes of channels K*(892)K+c.c. as 

M r C3f 7 o'" + = (g + g S )/2 + gf d , 
M^%- > K *- R+ = -e (l/2-s e ) + 3/2e 1 + ^, 

-o ( 15 ) 

ci + "- = (g+g s )/2+ g f, 

M*2%t cc - = -eo(l/2 - s e ) - 3/2 ei + a", 

as given in Tab. HI With those amplitudes we globally solve the Eq.fllJ) with all of the 
observed channels. Fixing sy= sp =0 and leaving all of the other parameters including g^ dd 
and free, the fit gives the most probable values of parameters |g|, g s , g^ dd , cos^o, |e |, 
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s e , a M |ei|, Si, r, dp. Those results are listed in the fourth column ("Solution 3") of Tab. 
Ulllwith x 2 / n do/=9.26/ll=0.84. From the parameters of solution 3 and Eqs.©,®,©> © 
and (1151) . we can get the production cross sections of the twelve VP channels, including the 
zero measurement of channel (pir . For example, for the channels pir and K*(892)K+c.c, 
the decay cross sections from ^(3770) resonance, and their decay branching fractions and 
partial widths can be calculated as 



a pn = 

" res 


(18.68itg) 


pb 


B( m = 


(1.87j££) x 


10~ 3 


ppvr _ 


51.1±ii;| keV 




K*~K 
V res 


= (0.56iJ;l) 


pb 


B K*K 


= (0.561JS) 


x 10" 




= (l-52l«g) 


keV. 



(16) 



ctm 



pv 

ctm 



And the production cross sections of the three pure E-M channels can be calculated as 

> 

pb 

(7.92 ±0.41) pb > (17) 
(2.69 ±0.29) pb. y 

The measured values in Eqs. (fT6l) . (|T7l) are all consistent with those measured in the last two 
subsections. As for the measurements of other channels, we have 



ctm 



B 



LOT] 



(2.1218:2) x io- 



BT> < 0.25 x 10 



-4 



at 90% c.l.) 



-4 



= (0.87±g;H) x 10 



B^' < 0.60 x 10~ 4 , ( at 90% c.l) 



(18) 



obtained in the assumption of that there exist the counter terms g^ dd and a M and sy=sp=0. 
If we fix the strong couplings g and g s to be the values which force the partial widths of pn 
and K*K etc. to be at the order of 1 keV which is at the same order of the J/ip VP decay 
coupling dynamics, the fitted x 2 is 47.27 for 15 degree of freedom which corresponds to 5.3 
standard deviation worse than that of solution 3. 

If there essentially were no the two counter terms specially for channel K*K+c.c (g^ dd = 
a M = 0), i.e. there existed the negative coupling g s and negative correction s e universally 
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allowed on all of other relative channels with the value s quark, the fit gives a somehow poor 
solution with large x 2= 19.2 for 13 degree of freedom, which gives 2.7 standard deviation 
away from the counter term assumption (solution 3). However, in this case, the effects of 
the opposite g s and negative s e in the related channels of (prj and (prf might be compensated 
by somehow larger nonet symmetry breaking with non-zero Sy and Sp. Ignoring the terms 
g^ dd and a M used in the last solution and leaving the nonet symmetry breaking parameters 
sy and sp free, we solve the Eq.(T4]) again to obtain a new solution. The results are listed in 
the fifth column ("Solution 4") of Tab. EH The fit gives x 2 /n d .o.f=13.2/ll = 1.2, which is 
a little bit higher than those in other solutions. However, the solution still has the statistical 
significance more than 5.6 standard deviation to the one obtained in the assumption that 
only the E-M components (i.e. the e , e x and d>i, and the rj — rj' mixing angle) act on 
the VP channel production. As for the solution 4, we obtain the unusual solution again, 
which is like solution 2 with the negative g s and negative E-M correction s e for the s quark 
couplings. The parameters of r, sy and sp which associate with the nonet symmetry breaking 
measurements are now with unexpected larger values. Nevertheless, as guessed above, those 
unusual numbers may suggest more complicate dynamics for the s quark production either in 
the channel K*K or in all of the relative channels, involving value s quark production. This 
solution predicts the decay cross sections, branching fractions and decay width of channels 
pir and K*(892)K and the production cross sections of three pure E-M production channels 
U7i°, prj and pr)', 



a pn = 

res 


(15.50^) pb 


B pn = 


(1.55^) x 10- 3 


Y p7T = 


42. 4^ keV 


n K*T< 

° res 


= (1.14iJi 2 7 ) pb 


B K*K 


= (1.14j£g) x 10 




= (S.lli 3 ;^) keV 


a ctm ~ 


= (11.671831) Pb 


n m - 

w ctm 


: (8.101SS) pb 


n PV' _ 
w ctm 


= (2.43°J 39 37 ) pb, 



(19) 



which are consistent with those obtained in the last solution with the counter terms. How- 
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and 07/', we obtain the branching fractions of 



(0.92^) x 1(T 4 



(20) 



ever, for the channels our], urj', <p\ 

B wr> 

B^' < 1.66 x 10" 4 , ( at 90% c.l.) 

< 1.57 x HT 4 , ( at 90% c.l.) 
B H = 3A6±lf 5 x 10" 4 , 

which are quite different comparing with those given in Eq. ffTSl) . As for the branching 
fractions for -0(3770) — > 07/, the measured values obtained by this analysis are also somehow 
inconsistent with that measured by CLEO-c Collaboration [6j, (see Tab. II Vj) . However both 
the results show the large '0(3770) VP decay couplings. The summed total cross sections 
over all VP channels of a^ 1 at both energy points of ^=3.670 GeV and 3.773 GeV are 
listed in Tab. IIVI The summed values of solution 3 and solution 4 as well as the results of 
CLEO-c measurement show in almost equals at the two energy points. The real resonance 
decays have been "hidden" . This is mainly owing to the destructive nature of interference 
between the couplings g and e . 



IV. DISCUSSION AND SUMMARY 

From above global analyses of the production cross sections for twelve channels of e + e~ — »■ 
VP measured by CLEO-c Collaboration at the -0(3770) resonance peak of yfs =3.773 GeV 
and at the energy of i/i=3.670 GeV in the continuum region, we obtained four different 
solutions of the parameters to control the VP channel production. From Tab. IIVI we 
see that the measured branching fractions or the production cross sections except for the 
channels involving the nonet symmetry breaking are all consistent with each other in the 
four solutions. 

It is remarkable that the large SU(3) symmetry strong decay strength |g| leads to the 
huge branching fraction of a level of 10 -3 for the typical channel pir, which corresponds to 
the decay width of two order of magnitude higher than that in J /if) decays. The large strong 
decay coupling hidden behind the VP channel production in the e + e~ annihilation at the 
0(3770) resonance peak might help people to understand the sources of the -0(3770) non-DD 
decays and give people some useful information to reexplain the long-standing pn puzzle in 
the 1 charmonium state VP decays. Furthermore, the large strong decay coupling with 
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TABLE IV: Summary of the results. The errors for our measurements are statistical only. The 
branching fraction of channel (prj measured by CLEO-c Collaboration is directly from work [6] and 



the upper 
of paper 6] 



imits of other channels measured by CLEO-c Collaboration are from the "Method II" 







Solution 1 


Solution 2 


Solution 3 


Solution 4 


CLEO-c 


pir 


o"3670 (pb) 


6.86^ 12 






6.43±i;|l 


8.0t 4 7 ± 0.9 




<y T (pb) 
a R (pb) 
Br 10~ 3 


5 35+ 1 ' 45 
, „, + 0.61 

i!iJ -o.67 




18.68l 4, 3g 
-, o 7 +0.46 
1 ' S7 -0.34 


5.08± 3 ; 2 7 

15 50+ 5 ' 04 
10 - ou — 2.96 
-1 -5+0.50 
i.DO_ 30 


4.4 ± 0.3 ± 0.5 

< 0.04 

< 0.004 


K'°k° + c.c. 


ct t (pb) 
a R (pb) 




W-30tf.ll 


19 34+ 1129 
la "° — 6.78 

0.28+g;S 


19.44 ± 3.61 

57+ 56 
u -° ' — 0.29 


23.3 ± 1.1 ± 3.1 
< 20.8 




Br 10- 4 




o.a»tS:5Z 


o.28±g;«t 


57+°' 56 
u -°'— 0.29 


< 20.8 




<y T (pb) 




< 0.502 


< 2.59 


< 0.77 


< 0.6 




art (pb) 




o.M±g;?S 


0.28+°;^ 


57+ 56 
u,D '-0.29 


< 0.1 




Br 10~ 4 




o. a8 ±g;!S 


n ,0 + 0.64 
0.28_ jj 


57+0-56 
u,D '-0.29 


< 0.1 




<y T (pb) 




11 74+0.31 
11. r^_ 34 


1 1 7Q + 0-33 
11. 'y_ .36 


11 67+ ' 36 


14.6 ± 0.6 ± 1.5 




a R (pb) 




0.0 


0.0 


0.0 


< 0.06 




Br 10~ 4 




0.0 


0.0 


0.0 


< 0.06 


A 


a T (pb) 




8.00 ± 0.22 


7.92 ± 0.41 


8.10 ± 0.44 


10.3 ± 0.5 ± 1.0 




trjj (pb) 




0.0 


0.0 


0.0 


< 1.3 




Br 10~ 4 




0.0 


0.0 


0.0 


< 1.3 


A' 


"T (pb) 




2.57 ± 0.28 


2.69 ± 0.29 


2.43 ± 0.32 


3.8±q 9 ± 0.6 




a R (pb) 




0.0 


0.0 


0.0 


< 0.4 




Br 10~ 4 




0.0 


0.0 


0.0 


< 0.4 


OJt] 


a T (pb) 
a R (pb) 


45+1-07 
2 19 +1 13 




n qo+0.62 
2.12 ± 0.21 


32+ 42 
u ""_0.10 

92+ 92 
u,y -0.77 


0.4 ± 0.2 ± 0.1 
< 0.1 




Br 10~ 4 


2.19 




2.12 ± 0.21 


92+ ' 92 
u - a — 0.77 


< 0.1 


U1T)' 


<y T (pb) 


44+ - 46 




59+ 0B9 

U.OM_ Q 3g 


44+ 1 - 84 
u, *^-0.41 


0.6ig'g ± 0.6 




"R (pb) 


< 0.33 




< 0.27 


< 0.47 


< 1.9 




Br 10~ 4 


< 0.33 




< 0.27 


< 0.47 


< 1.9 


4>V 


"T (pb) 


3.99lj 27 




3 91+ 275 


3 70+ 2 7s 
J ' " J -1.02 


4.5 ± 0.5 ± 0.5 




o- H (pb) 


n 84+ ' 79 






< 1.45 


2.4 ± 0.6 




Br 10- 4 


n 84+ ' 79 




o.87i8;SS 


< 1.45 


3.1 ± 0.7 


4>v' 


a T (pb) 


1.95_ 2 g3 




2 29+ 4 20 


1 69+ 1 ' 37 


2.5±J l ± 0.4 




o"ij (pb) 


< 0.52 




< 0.66 


3 46+ 4 ' 90 
J '* D -1.76 


< 3.8 




Br 10~ 4 


< 0.52 




< 0.66 


3 46+ 4 90 
J '* D -1.76 


< 3.8 


All "eft," summed 


CT 3670 (P b ) 






58.3 ± 6.2 


62.3 ± 6.9 


64 o 4 " 158 
-6.8 




°T,3770 (pb) 






54.4t 7 3 4 ° 


52.9+^; 3 


64.61 4 ;! 



the opposite s-quark strong coupling and the minus E-M correction s e or equivalently, the 
s-quark "anomalous magnetic moments" a M , required by the K*K production as presented 
in those solutions with different treatments and different assumptions in this work are all 
unusual comparing with the conventional hard gluon annihilation picture plus single ?/>(3770) 
resonance assumption. If those measurements and analyses are all correct, one has to re- 
understand the strong interaction dynamics which leads to the large Okubo-Zweig-Iizuka 
rule violation in the vector meson -0(3770) decays and the strange behavior of the s quark 
couplings to the light hadron production in the energy region around ^(3770) resonance. It 
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seems that people has to seriously consider the role of the long distance strong interaction 
corrections including the D(D S ) meson exchange scheme to describe those anomalous phe- 
nomena and the large Okubo-Zweig-Iizuka breaking in this energy region, like many authors 



did 



is, 
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The destructive interferences mechanism might manifest the possibility of the significant 
existences of buried ^(3770) non-DD decays. However, how do people understand the 
large net inclusive non- DD hadron branching fraction of ^(3770) decays measured by 
BES Collaboration {3, 4] recently. Phenomenologically, for example, if the behavior of the 
reversed s quark strong coupling appeared in VP channel is still maintained in those cases, 
it would lead to constructive interference with the parallel continuum E-M amplitudes and 
cause abundant strangeness meson production at the ?/>(3770) peak, resulting in the large net 
cross section excess. The one of the exceptions is the channel K* K+c.c, in which the E-M 
production is due to the magnetic moments coupling. The minus E-M coupling correction 
for s quark or the counter term a M as a special "anomalous magnetic moments" enhances 
the E-M production at the continuum region and leads to the observation of the equal cross 
sections of channel K* K+c.c. at the two energy points -\/i=3.670 and 3.773 GeV in work 
js]. This argument can be cleared up by coming more precisely experimental measurements. 
Of course, another probable outlet would be that there are more complicated structures 
or contents in the ^(3770) resonance scope which are evident in a measurement of cross 



sections for e + e~ — > hadrons by the BES Collaboration [20|, |2l|. This measurement of the 
cross section indicates that there are somehow complicate "diresonance" structure instead 
of the conventional single ^(3770) resonance assumption 20|, |22|. The exist of the extra 



substances might respond to the unusual behavior ofj/> 
large non-DD branching fraction of ^(3770) decays [3j, 



3770) VP decays and the measured 
4|- 
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